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Direct Illumination 
Rendering equation: Lout -| LBVdw 
[Kajiya86] Q 


L= Lighting e 


Monte Carlo: 


Lout = FL LBV 


Lout 


N B= BRDF x cosine 


Importance Sampling 


Send rays in directions with large LBV 
But, V is unknown so hard to include in sampling 


occluder 


LB (Light x BRDF) sampling Our algorithm 


Related Work 


Kollig & 


Debevec Cohen & 
Keller ostromoukhov 


[1998] Debevec 


Veach & 


Ward Cani n 
Donikian 
рд ршн Ashikhmin Caustres бий onikan быша (Stab 
BRD = talonde [2000] (озу 12006] 
i e 
a rn cad 1 
йт «tal. БЛ 120061 
E 
pou 
Tee 
man ^ш 8002) 
Er 
11538) 


Time 


Outline 


1. Visibility Measurements 


Statistical correlation of the visibility function 


2. The Visibility Cache 
Similar to radiance caching, but for visibility 


Applications: 


= Reduce noise in Monte Carlo ray tracing 


= Ambient occlusion 
= Fast lighting preview (aka PRT) 


1. Measurements 


Large number of random pairs of positions 


d = Distance 
@ = Normal difference 


Visibility Correlation 
We see the visibility maps as distributions 


Compute the correlation coefficient: PX,Y 
(Pearson’s product-moment coeff.) 


Examples: (uà "Ч -1 


Scene (a) 


Correlation 


Norma! difference (deg) 


Correlation 


Scene (a) 


ин 
Distance 


Scene (c) 


Correlation 


Normal difference (deg) 


as 


90 


135 


Correlation 


Scene (с) 


15 2 28 3 38 
Distance 


Scene (d) 
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2. The Visibility Cache 


During rendering, insert cache records: e 


Each record is a 
low-res visibility 
map (e.g. 32x32) 


Typically scene: 
10k - 20k records 


= 5-10 M rays 


Where to place records? 


Several heuristics based on measurements 
a G j 
eiai Details in 
= Average visibility difference the paper 


= Max number of wrong visibility queries 


Example: Static Scene 


Example: Static Scene 


Video 


| Example: Dynamic Scene 


Example: Dynamic Scene 


Implementation 


00000111010111111010110110...10110000000 


Implementation 


0000011101011111 


1010110110...10110000000 


Implementation 


000001110101[111111010110110...10110000000 


Visibility Difference 


0000011101011111101011011011001011111111111000001111110110000000 
0000001101111111111000001011000011101101001001001111111110100100 


00000100001 0000001 0011010000001 0000100101 10001000000001000100100 


1) XOR 2) Count bits = 15 = 23.4% difference 
SSE: 128 bits at the time also using SSE 


Memory Use 


Example 


Visibility map size: 32x32 pixels 
1 bit/pixel > 128 bytes + 44 bytes book-keeping 


Number of cache records: 11,900 
Total memory use: 2 MB (fits in CPU’s L2 cache) 


App #1 Monte Carlo 


Use Control Variates to reduce noise: 
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Original Correlated Difference has 
function approximation lower variance 


Da / LBVdw = / LBV - LBVdw + / LBV dw 


App #1 Monte Carlo 


Use Control Variates to reduce noise: 
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function approximation lower variance 
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Given by Visibility Cache 


App #1 Monte Carlo 


Use Control Variates to reduce noise: 
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Given by existing product sampling algorithms 
[Cline et al. 2006, Clarberg et al. 2008] 


App #1 Monte Carlo 


Use Control Variates to reduce noise: 


fe кх) fed-se) 
I 
Original Correlated Difference has 
function approximation lower variance 


Unbiased! 
Fe / LBV dw = / LBV — LBVd + / LBVdw 


Evaluated using Given by quick 
importance sampling quadtree product, 


Equal-Time Results 


State-of-the-art Our algorithm 
product sampling (30 samples/pixel) 
(39 samples/pixel) 


Product Sampling 


Our Algorithm 


Product Sampling 


Our Algorithm 


Product Sampling 


Our Algorithm 


App #1 Summary 


Good use of control variates (finally) 
Unbiased even if visibility approximation is off 
Animations: Re-use over time without lagging 


Up to 5x variance reduction 
Depends on scene complexity 


= Some scenes have too low visibility correlation 


= Largest improvement is in shadow regions 


App #2 Ambient Occlusion 


Pre-integrate records: 4; -f Vi )m- )d 
a 
Approximate as weighted sum: A(x) = YO w;A; 


Issues 


Amortized cost per pixel is very low, 
but the quality is not (yet) perfect 


Hard to detect small geometric features: 


2 n 4 2 Same here 


Big risk we miss the step when N shadowing 
sampling the hemisphere 


Issues 


Amortized cost per pixel is very low, 
but the quality is not (yet) perfect 


Hard to detect small geometric features: 


= We use minimum hit distance [Tabellion & Lamorlette 2004] 


5510 lo 1 


Exaggerated blurring due to filtering, some blotchiness 


= We expect “visibility gradients” (~ irradiance gradients) to help 


App #3 Fast Lighting Preview 


Similar to Precomputed Radiance Transfer 
Cost: One sparse quadtree traversal / pixel 


Conclusion 


Visibility is often strongly correlated 
But, some scenes are just too complex (e.g. grass) 


Visibility Caching 
Store compact visibility maps sparsely in the scene 
Interpolate and use for: 

= Monte Carlo rendering using control variates 


= Ambient occlusion 


= Lighting design 
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